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The primary goals of this thesis included developing a computer-integrated LDV 
system, utilizing the LDV system to measure the mean and fluctuating quantities of the 
flow field of an airfoil-spoiler system, and comparing the results of the experiment with 
previous hot-wire anemometry results. 

The experiment was conducted in the NASA Ames 1 1 cm x 25 cm indraft 
tunnel. A two-component LDV system was developed and set up, and measurements 
were obtained along the upper surface of the airfoil (ahead of and behind the spoiler) 
as well as in its wake, close to the trailing edge. Mean velocity profiles and turbulent 
stresses were calculated and comparisons were made with previous hot-wire 
experiments. 

It was shown that the hot-wire technique is not as accurate as the LDV in highly 
turbulent regions of the flow field. The hot-wire results are also not reliable in the near 
wake region due to the inherent reverse flow. In other, less turbulent, regions of the 
How field where there is no reverse flow, good agreement between previous hot-wire 
data and present experimental results was found. 
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I. INTRODUCTION 



A. SPOILERS 

A spoiler is a plate on the upper surface of a wing which can be deflected to spoil 
the flow and to cause early separation. Consequently, lift decreases and drag increases. 
Spoilers are used as effective aerodynamic control surfaces in transport aircrafts. They 
can be symmetrically deflected as speed brakes or as lift dumpers at touch down. 
When spoilers are deflected asymmetrically, they can provide substantial amounts of 
roll control in certain conditions of flight. 

The flow field of an airfoil with a deflected spoiler is complex: including flow 
separation, reattachment and vortex shedding [Ref. l:p. 1]. 

The basic steady state structure of the flow field of an airfoil with a deflected 
spoiler has been measured and hypothesized by the Boeing Research Group [Ref. 2]. 
Further work into how the unsteady wake relates to the airfoil geometry and to other 
parameters was undertaken by Lee, ct al. In their experiments, Lee, et al., used 
hot-wire anemometers and pitot tubes to systematically measure the mean and 
fluctuating quantities on the airfoil surface, as well as in the wake, in order to better 
understand the structure of this flow field. [Ref. 1] 

B. LASER DOPPLER VELOCIMETER SYSTEMS 

Advanced research in experimental fluid dynamics requires a familiarity with 
sophisticated measurement techniques. Optical methods and in particular, the Laser 
Doppler Velocimeter (LDV) are now recognized as the most reliable means for 
performing measurements in complex turbulent flows [Ref. 3]. 

Since its inception in 1964, the LDV technique has been revised, developed, and 
tested by many researchers and its capabilities are reasonably well-understood. The 
method has the advantage of performing instantaneous, non-intrusive velocity 
measurements of small particles suspended in the flow. These measurements are 
relatively independent of fluid properties and the method does not require frequent 
calibration. Furthermore, its ability to measure the flow without disturbing it is 
especially attractive. However, because of the range of disciplines involved including 
optics, electronics, light scattering, and signal processing, the implementation of this 
method can be quite onerous. Serious measurement errors can occur if the 
measurement techniques are not applied correctly. 
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C. THESIS GOALS 

Considering this previous work, the goals of this thesis experiment were threefold. 
First, to develop and set up a two-component LDV system with complete 
computer-experiment integration employed. Second, to utilize this LDV system to 
measure the mean and fluctuating velocities of an airfoil with a deflected spoiler flap, 
especially in the reverse flow region and the near wake where the fluctuations are large. 
Third, to compare the LDV data obtained with the previous hot-wire anemometry 
data. 
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II. EXPERIMENTAL ASPECTS 



A. EXPERIMENTAL SETUP 

The experiments were conducted in the NASA Ames 1 1 cm x 25 cm indraft 
wind tunnel. Figure 2.1 shows the test section of this tunnel. The test section is 
equipped with slotted walls and a number of plenum chambers above and below the 
centerline. These features were used in previous experiments to give the tunnel 
adaptive walls. The adaptive wall concept was not employed in this experiment, 
however the slotted walls did provide a means of reducing wall interference. The wind 
tunnel also is equipped with interchangeable glass and plexiglas test section windows. 
For these experiments, the glass windows were used due to their higher refractive 
qualities which enhance the LDV signals. 

The basic model used was a Boeing advanced transport research airfoil section 
with a 6-inch chord and maximum thickness ratio of 11.3% as shown in Figure 2.2. 
This airfoil has a spoiler of 15.5% chord length hinged at 73% chord. The coordinates 
of the airfoil are given in Table I. This airfoil was selected because of the availability 
of hot-wire anemometer and pitot tube experimental data for comparison with the 
two-component LDV results. 

B. INSTRUMENTATION 

1. Pressure Instrumentation 

For this experiment, pressure information was obtained by taking advantage 
of a bank of six, 2--port Scanivalves tapped into the tunnel during a previous 
experiment. The six pressure transducers were of the differential type with a range of 
±5 psig. One side of each transducer was connected with the measured pressure 
through plastic tubing while the other side of each transducer was vented to 
atmosphere. 

In order to calibrate the transducers, one port of each Scanivalve was tapped 
to a calibration pressure generated by a small electric vacuum pump. Another port of 
each Scanivalve was tapped to the atmosphere. The calibration pressure and 
atmospheric pressure were each measured by alternate methods. These two known 
pressures and the values given for them by each transducer allowed the Scanivalves to 
be calibrated. 
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Figure 2.1 NASA Ames indraft wind tunnel test section. 
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Figure 2.2 Boeing advanced transport research airfoil section. 

2. Laser Doppler Velocimeter Instrumentation 

Figure 2.3 shows the two-component, two color LDV arrangement used at 
NASA Ames. The laser is a four Watt Argon ion unit (Spectra Physics Model 64). 
This laser is configured with an end cavity multiline reflector which allows it to operate 
in the multiline mode. Translation stages in three axes permit positioning of the probe 
volume anywhere in the test section. 

In the multiline mode, the resultant beam is comprised of all wavelengths 
inherent to the Argon ion unit. This beam exits the laser in a fairly collimated state, 
hence an external collimator is not required. The beam first passes through a 
half-wave plate (polarization rotator) where it becomes horizontally polarized. A set of 
Brewster angle prisms is then used to separate the primary beam into its various 
wavelengths. To obtain maximum light transmission through the prisms, the beam 
polarization must be parallel to the plane defined by the incident and refracted beams 
[Ref. 4]. This was the reason that the initially vertically-polarized laser beam was 
passed through a polarization rotator prior to entering the dispersion prism. When the 
two prism method is used, the prisms also turn the beams which, with a little planning, 
may be used to advantage. At some distance from the prisms, the beams have 



15 



TABLE I 

AIRFOIL COORDINATES (INCHES) 


X 


Y 

upper 


y 

lower 


0.0 


0.0 


0.0 


.0030 


.0222 


-.0108 


.0060 


.0300 


-.0162 


.0150 


.046S 


-.0258 


.0300 


.0672 


-.0348 


.0900 


.1242 


-.0588 


.1200 


.1452 


-.0672 


. 1 500 


. 1626 


-.0750 


.2250 


.2004 


-.0912 


.3000 


.2304 


-.1050 


.4500 


.2790 


-.1296 


.6000 


.3126 


-.1524 


.7500 


.3396 


-.1728 


.9000 


.3594 


-.1920 


1.2000 


.3858 


-.2250 


1.5000 


.3990 


-.2502 


1.S000 


.4044 


-.2670 


2.1000 


.4038 


-.2748 


2.4000 


.4002 


-.2742 


2.7000 


.3924 


-.2658 


3.0000 


.3816 


-.2502 


3.3000 


.3678 


-.2298 


3.6000 


.3510 


-.2064 


3.9000 


.3288 


-.1818 


4.2000 


.3018 


-.1560 


4.5000 


.2700 


-.1308 


4.S000 


.2304 


-.1044 


5.1000 


.1842 


-.0792 


5.4000 


.1308 


-.0540 


5.7000 


.0696 


-.0282 


6.0000 


.0024 


-.0024 



sufficient separation so mirrors can be used to direct the individual beams to the next 
optical components. The two stronger lines, green (514.5 nm) and blue (4S8 nm), are 
guided by means of turning mirrors to Bragg cells. The weaker, unused lines are 
"trapped" (blocked by nonreflective plates) to keep them from interfering with the 
system. Prior to reaching its Bragg cell, the green beam is passed through another 
half-wave plate in order to vertically polarize it. 

A Bragg cell is an acousto-optic device. Both the blue and green beams are 
passed through Bragg cells which split each into two equal-intensity, slightly diverging 
beams while inducing a 40 MHz modulation in one of the beams in each pair. The 
difference is that the green beam is split in the vertical plane while the blue beam is 
split in the horizontal plane (Figure 2.3). 
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Figure 2.3 Two-component LDV arrangement used at NASA Ames. 
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In each case (blue and green), the two diverging beams are then made parallel 
by passing each pair through its own alignment prism. From there on, turning mirrors 
guide the beams through the translation stages. They are finally focused by means of a 
condensing lens at the desired point in the test section of the wind tunnel. 

The four beams must all focus to a common point in the How and be 
approximately coaxial. Using the steering mirrors for the appropriate beam pair, the 
focused beams are adjusted to visually overlap. Following this, various techniques 
exist for making fine adjustments to insure the beams all cross at a common point. 

An off-axis, forward scatter configuration has been selected for the receiving 
optics used at NASA Ames. This gives flexibility to the system so it can be adjusted to 
different requirements; if maximum signal strength is required, angular deviation from 
the beam axis of symmetry is minimized. If focusing on the smallest part of the probe 
volume is of interest, the angular deviation is made as large as the physical constraints 
will allow. 

The first lens in the receiving optics is the aperture through which the 
scattered light is collected (Figure 2.4). The second lens collimates the beam. A 
collimated beam is required by the polarization beamsplitter for maximum efficiency. 
The polarization beamsplitter allows the vertically polarized green beam to pass 
straight through while the horizontally polarized blue beam is reflected upward. The 
use of polarization rather than frequency to separate the two wavelengths substantially 
increases the efficiency of the system. Fach beam then passes through a line filter 
centered on the respective wavelength (blue or green) in order to reduce both the 
ambient light and the cross-talk due to impure polarization. Each line filter is followed 
by a lens which serves to focus the light onto the pinhole aperture of the respective 
photomultiplier. The photomultipliers convert this light energy into electronic signals 
which are then passed on to the signal processors. 

Preamplifiers are located in the photodetector housings to further increase the 
signal amplitude and drive the signal over relatively long cables. The preamplifier is 
kept close to the detectors in order to increase the signal level before additional noise 
enters the signal cables. It is always good practice to amplify the signal as much as 
possible with the photomultiplier tube before going to the preamplifier which inevitably 
produces electronic noise. 

The signal processors used in this experiment were of the counter type and 
produced by Macrodyne. These processors have good resolution and superior noise 
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Figure 2.4 Computer-experiment integration. 
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rejection circuitry. The Macrodyne unit is a high-speed counter with a 1 GHz clock 
designed to operate in a variety of experimental situations. It has been applied under 
single particle, multi-particle and multi-velocity conditions. In addition to the basic 5/8 
(or 10/16) signal validation criteria and oversize particle rejection, a sophisticated 
multi-level sequence validation has been incorporated. This allows the experimenter to 
select the criteria which will best yield data for a given signal-to-noise ratio, signal 
level, seeding level, Doppler frequency, and fluctuation, indicative of the test situation. 

When the particle seeding rate is low, as w r as the case in this experiment, the 
LDV signals become intermittent bursts. Under these conditions, the counter 
processor is the best choice as a means for obtaining the Doppler difference frequency. 
If the seeding is high, the signals will be almost continuous with random amplitude and 
phase fluctuations. The counter will still function but cannot be expected to produce 
accurate results to better than ± 1/2 cycle. 

The counter processor determines the Doppler difference frequency in the time 
domain. Timing over a selected number of cycles in the burst signal is accomplished 
with a crystal controlled "clock" oscillator. The frequency of the oscillator must be 
much higher than the Doppler difference frequency to achieve good resolution. The 
clock count is initiated when the high-pass filtered signal exceeds a set threshold level. 
A preselected number of cycles (8 or 16) are processed and used to stop the counter. 
Greater accuracy may be obtained by processing more cycles. Some processors, (e.g., 
Aerometrics) count all of the cycles in the burst and average the result over the number 
of cycles counted. This capability is especially important when the seed particles are 
polydisperse and produce a range of signal amplitudes. In such cases, there is the 
possibility of completing the preset number of cycles on the early part of the signal 
which has a lower signal-to-noise ratio and slew rate. There is also the possibility of 
reading signals from larger particles more than once. 

Although the importance of attaining good quality signals cannot be 
overemphasized, in many practical cases the signals will be noisy. Therefore, each 
burst must pass certain criteria before it is accepted. The so-called periodicity check is 
used to reject: 1) signals with insufficient cycles, 2) signals produced by two or more 
particles in the measuring volume at a time, and 3) noise. Two counters are used to 
make period measurements over a different number of cycles, usually 5 and 8. The 
ratio of the two measurements is then compared to five-eighths and if they agree to a 
specified percentage, the signal is validated. The problem with this technique is that 
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the entire signal may be rejected because of a noisy cycle early in the burst. However, 
when the method is coupled with the three-level validation circuitry, it can be effective 
in rejecting spurious noise. 

The three-level validation requires that a signal. Figure 2.5. exceed a positive 
threshold level, pass through zero, and exceed a negative threshold level before the next 
zero crossing is counted as a cycle. Such logic has proven to be very effective in 
rejecting noise. If the conditions are not met on any cycle, the counter resets and 
starts again almost immediately. Thus, noisy cycles on the start of a burst need not 
cause the loss of an otherwise good signal. 

These amplified electronic signals are sent to the signal processor and 
ultimately to the data processor for conversion into velocities. 

Ideally, an LDV needs no calibration since it provides a direct measure of the 
particle velocity. The fringe spacing. 6, can be determined using the following formula 
where /. is the laser wavelength and 0 is the angle of the intersecting beams: 



6 = /. (2 x sin(0 2)) 



Knowing the fringe spacing and the time required for a given number of crossings, 
panicle velocity can be determined. In practice this method can be inaccurate due to 
misalignment of the optics. The beam intersection angle, 0, is also very difficult to 
compute accurately. Therefore, it is wise to test the entire system by measuring the 
light scattered from an object moving at a velocity that is known a priori. A rotating 
disk with a mark at a known radius is usually the most convenient means of doing this. 
However, another method (the method used in this experiment) is to measure the free 
stream velocity using the LDV and compare this value with that computed from the 
Scani valve data. 

3. Flow Seeding Requirements 

An inconvenient characteristic of the laser Doppler velocimeter technique is 
the need for seed particles. The LDV actually measures the velocity of small particles 
that are assumed to move with the flow. If these particles do not have the appropriate 
size and concentration distribution, the experiment results may be adversely affected. 
Particle sizes on the order of 1 micron are often quoted as adequate for most gas Hows 
and can respond to turbulence frequencies exceeding one kHz. Smaller particles may 
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be below the detectable limit of the instrument but will contribute to the background 
noise while larger particles that are easily detectable may not adequately track the flow. 

It was determined in this experimental situation that seeding of the flow was 
not absolutely necessary. Several runs were made with no seeding supplied. The data 
obtained in these runs was accurate as dust particles acted as the flow-following agent. 
However, it was determined that by using a flow seeding system the data rate was 
increased by a factor of four to five. The system used consisted of a paint sprayer 
spraying latex particles (on the order of 3 microns) suspended in alcohol. 

C. COMPUTER-EXPERIMENT INTEGRATION 

The experiment was fully integrated with a Data General NOVA minicomputer 
for on-line data acquisition and analysis, Figure 2.4. A main program 
("SPOILER. FR") was written to. perform all tasks and record all data required during 
the experiment so that individual tunnel runs proceeded virtually 'hands off'. The main 
program called a host of other programs during each run in order to accomplish 
various tasks. These other routines are listed here: 

(1) "SCANT. FR" - used to obtain pressure data from the Sea ni valves and 
compute free stream velocity and free stream Mach number. "SCANT. FR" 
itself calls three other routines in order to determine ambient temperature, 
calibration pressure and barometric pressure: 

(a) "TEMP.FR" - which vields ambient temperature from a thermocouple 
via an analog-to-digital converter. The thermocouple is located 
upstream of the test section. 

(b) "CALPRES.FR" - which calculates the output pressure from the 
calibration vacuum pump. The pressure value is read by this 
program using an analog-to-digital converter. This value is then used 
in the calibrafion of the Tscanivalves. 

(c) "NPARO.FR" - which yields the value for current atmospheric 
pressure from a Paroscienfific pressure sensor via a digital I/O. 

(2) "STDTM.FR" - which gives the current time and date in order that the data 
files being created by the main program can be labeled. 

(3) "TRELAY.FR" - to command the taskmaster to move the traverse mechanism 
to the next data collection location. 

(4) "PLOTP.FR" - to provide on-line data analysis by allowing any of five curves 
to be plotted immediately after a given run. ‘ 

D. EXPERIMENTAL CONDITIONS 

The free stream velocity was maintained at a nominal value of 70 m/sec, 
corresponding to a Reynolds number of 7.4 x 10 5 based on the airfoil chord. 
Experiments were carried out with the spoiler deflected 30 degrees and the angle of 
attack at 0 and 5 degrees. Locations where data was obtained are shown in Figure 2.6. 
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Figure 2.6 Data collection locations. 
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In order to determine the velocity profiles, data was acquired at varying steps 
based on the velocity gradient. Fine step sizes were taken near the airfoil for boundary 
layer profiles and near free shear layers in the wake measurements. Typically, 15 data 
points were taken in the boundary layer and 20 points in the wake. The 4-beam LDV 
arrangement precluded obtaining data closer than 0.1 inch from the airfoil surface, 
however. 

The boundary layer traverse was along the vertical axis of the test section, so the 
traverse direction is not normal to the airfoil surface as are usual boundary layer 
profiles. Angle of attack and model surface curvature make the measurement axis 
incline to the model. 
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III. DATA ACQUISITION 



A. FREE STREAM VELOCITY 

At the beginning of every run, the main experiment-driving program 
("SPOILER. FR") calls for the Scanivalve pressures, ambient temperature, calibration 
pressure and atmospheric pressure to be read. Assuming isentropic flow in the test 
section upstream of the airfoil and away from any shock waves, the free stream Mach 
number, M, can be calculated from the stagnation pressure, P 0 , and the static 
pressure, P, along the wall: 



M = {((Pq/p/Y-O/Y- 1)(2/(y-D)} 1/2 



For y = 1.4: 



M = {((P 0 /P)^ 2 / 7 > — 1) x 5} 



Then, velocity is given by: 



V = Ma 



where 



a = VyRT 



T is obtained from the total temperature, T Q , using the following formula: 



T = T 0 /{1 + ((y-l)/2) x M 2 } 
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T q is measured by a thermocouple at the wind tunnel inlet. 

The above method is used to obtain free stream Mach number and free stream 
velocity for every run. The value for free stream velocity is used in the data reduction 
process to normalize the airfoil velocities obtained from the LDV. 

B. AIRFOIL VELOCITIES 

Processing of the Laser Doppler Velocimeter data obtained during this 
experiment was carried out using the same Data General Nova digital computer and 
driving program used for all other experiment functions. "SPOILER. FR", the main 
program for this experiment, calls a digital I/O driver program ("DI024") to read 
words as fast as possible from two cards (these cards receive the digital output from 
the two Macrodyne channels, one for each LDV component). As soon as 1000 
samples (per component) have been taken by the Macrodyne processor and read by 
DI024, the binary output from each sample is converted to decimal using a decoding 
subroutine ("DECODE!. FR"). 

As mentioned earlier, the Macrodyne signal processor used in this experiment 
utilizes a 1 GHz counter. The number of "counts" required for the passage of 8 cycles 
is sent to the digital I/O in the form of a 10-bit mantissa and 4-bit exponent, namely 



D q , D s ,...D n x 2 ex P 



For example, if 400 "counts" were required for a particular Doppler burst, the 
output would be 



0000110010 0011 
Mantissa EXP 



which, upon decoding, reduces to 

50 x 2 3 = 400 
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Since the clock counts are produced at 4 GHz (the 1 GHz clock is divided into 4 
phases), each clock count is 0.25 nanoseconds. Thus, 



t 



s 



D 9 , D gf ...D 0 x 2 ex P x 



(025) 

( 8 . 0 ) 



ns 



is the period of the signal when 8 cycles are measured. 

Converting this to signal frequency, f s , yields 

1 2 5 ' cx P x 1000 

f = — = MHz 

b t $ Mantissa 

Since the signal frequency, f s , is composed of the Doppler frequency, the shift 
frequencyand the mixing frequency, the Doppler frequency can be obtained from 

^D “ — ^shift + ^rnix 

Once the Doppler frequency of a particle is known, its velocity can be determined 

from 



velocity = const x f^ 



where "const" is determined based on the fringe spacing or direct calibration. 

With the 1000 samples from each component converted to decimal, the data can 
be reduced to average and rms velocities using the formulas found in Table II. This 
data (both raw and reduced) is then stored for future analysis. 

While most data analysis is performed after completion of the experiment, several 
on-line analysis routines are available for use as real time performance checks. For 
example, "HISTOGRAM. FR" is a program available to first read a given number of 
samples from one of the Macrodyne channels and plot frequency versus the number of 
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samples at each frequency in histogram fashion. With this program it is easy to 
determine whether or not the data is conforming to the Normal Distribution (to which 
ideal data should conform). Another on-line analysis tool which is built into the main 
program "SPOILER. FR" is the subprogram "PLOTP.FR". This routine gives the 
experimenter the option to have any of five velocity profiles (i.e. , U^g/Ur*-,, 
V ave /Voo. etc.) plotted on the screen immediately after completing an experimental 
run. This subprogram was used extensively to immediately inspect the quality of each 
run. 



TABLE II 

VELOCITY PROFILE AND TURBULENT STRESS EQUATIONS 



1 N 

V= — Tu- 



»=i 



V = — Tv- 

*k" 



1 i = 1 



1 N 

v 2 = — y v 2 - v J - 

n ^ • 1 



i=i 






N t=i 
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C. REYNOLDS SHEAR STRESSES AND TURBULENT STRESSES 



Table II shows the equations used to determine mean velocities, turbulent 
stresses and Reynolds shear stress throughout this experiment. The value for X was 
generally 1000 since this was the number of velocity samples taken at each 
measurement location (for each component). 
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IV. RESULTS AND DISCUSSION 



A. INTRODUCTION 

One of the purposes of this experiment was to demonstrate the benefits of the 
two-component Laser Doppler Velocimeter by employing this type of system to 
measure the flow field of an airfoil with a deflected spoiler. In order to present the 
results of these measurements effectively, a graphics program called "PLOT-3D" was 
used to plot the data obtained. "PLOT-3D" is available on the VAX system at the 
NASA Ames Research Center. 

A data interpolation routine written by Chyang Sheng Lee of Stanford University 
was also used to interpolate the measured data into a grid system so that it could be 
plotted. All of the contour plots presented use the interpolation routine to interpolate 
in both the x and y directions. However, some of the velocity vector plots use 
interpolation in the y direction only. This was done because, in some cases, linear 
interpolation of the data gives a false impression of what the How field is actually 
doing. 

In all cases, the data presented in this section is for a = 0 s and 5 = 30 5 . 

B. FLOW FIELD ANALYSIS 

Velocity measurements were taken at a variety of points along the upper surface 
ahead of and behind the spoiler. Figure 4.1 shows a plot of the data obtained over the 
rear portion of the airfoil. Note the low-velocity, recirculating flow region directly 
behind the spoiler. 

While the data for Figure 4.1 was interpolated in the y direction only, the data 
for Figure 4.2 was interpolated in both the x and y directions. It can be seen in Figure 
4.2 that interpolation in the x direction yields some fairly unrealistic velocity vectors 
due to the relatively low number of measurement locations as compared to the number 
of grid points. For example, those directly ahead of the spoiler. On the other hand, 
complete interpolation (in both the x and y directions) does serve to paint a good 
overall picture of what the flow field is doing based on a finite number of velocity 
measurements. 
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Figure 4.1 Flow field measured over rear portion of airfoil 
(interpolation in y direction only). 
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Figure 4.2 Flow field measured over rear portion of airfoil 
(interpolation in both directions). 
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C. WAKE 

The wake behind the spoiler was measured with the two-component LDV system 
at numerous locations to determine the mean and rms profiles. The results from the 
wake arc presented in this section. 

1. Wake Profiles 

The mean longitudinal and transverse wake profiles measured by the LDV for 
a = 0° and 5 = 30° are presented in Figures 4.3 and 4.4. The symbols are the 
measured data and the lines represent the best fit curves of the data. As expected, the 
width of the wake, the velocity defect and transverse velocity all decrease in the 
downstream direction. 

The Reynolds shear stress, longitudinal and transverse turbulent stresses 
measured by the LDV are presented in Figures 4.5 to 4.7. Obvious in all three of these 
figures is the fact that the fluctuating region appears to grow from the x/C= 1.0 
(trailing edge) location to x/C= 1.25 then decrease as it approaches x/C= 1.5. The "S" 
shape profile in the Reynolds shear stress is typical for a two-dimensional wake 
[Ref. l:p. 29]. The longitudinal and transverse turbulent stresses show two peaks at the 
locations where the two free shear layers (corresponding to the airfoil and the spoiler) 
oscillate in the near wake. 

Figures 4.8 and 4.9 show the velocity vectors in the wake. A scaled profile of 
the airfoil is superimposed in each figure for reference. The data in Figure 4.8 w r as 
interpolated in the y direction only. Figure 4.9 was interpolated in both directions and 
depicts the mean wake velocity profile quite well. 

2. Contour Plots 

Figures 4.10 to 4.16 depict "contour" plots of various quantities in the airfoil 
wake. Figure 4.10 shows regions of equivalent U-magnitude (U is referenced to free 
stream velocity). Note that there are regions where the ratio U/Uqo is greater than 
unity. This is due to tunnel blockage effects which cause the flow to accelerate more 
than normal around the airfoil. 

Figure 4.11 depicts the regions of equivalent V-magnitude (also referenced to 
free stream velocity). A minus sign on a contour line in Figure 4.11 implies that the V 
component of velocity is negative (downward) in that region. 

In Figure 4.12, the total velocity magnitudes in the wake are contoured. 
Again, contour lines with values greater than one arise due to blockage effects. 



34 



x/C=1.0 x/C=1.25 x/C=1.5 



O 





n w o 

o o o o o 

I i 



oA 

Figure 4.3 Longitudinal velocity profiles in the wake. 
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Figure 4.4 


Transverse velocity profiles in the wake. 
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Figure 4.5 Reynolds shear stress in the wake. 
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Figure 4.6 Longitudinal turbulent stress in the wake. 
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Figure 4.7 Transverse turbulent stress in the wake. 
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Figure 4 9 Velocity vectors in the wake 
(interpolation in both directions). 
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Figure 4.11 Contour plot of V-magnitude in the wake. 
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Figure 4.12 Contour plot of total velocity magnitudes 
in the wake. 
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Figure 4. 1 3 Contour plot of flow angles in the wake. 
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Figure 4.14 



Contour plot of longitudinal turbulent stresses 
in the waRe. 
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Figure 4.16 Contour plot of Reynolds shear stresses 
in the wake. 
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Figure 4.13 shows how rapidly the flow angles change in the near wake (the 
closer the contour lines, the steeper the gradient). 

Longitudinal and transverse turbulence and Reynolds shear stress are 
contoured in Figures 4.14 to 4.16. Figure 4.14 clearly shows the two regions where the 
turbulence intensity is the highest. These two regions are probably caused by 
counter-rotating vortices originating at the spoiler trailing edge (clockwise vortex) and 
the airfoil trailing edge (counterclockwise vortex). 

D. COMPARISON WITH HOT-WIRE ANEMOMETER DATA 

One purpose of this thesis was to compare the data obtained in this experiment 
with the hot-wire data obtained under similar conditions by Lee [Ref. 5]. Figures 4.17 
to 4.21 are LDV/hot-wire comparisons of longitudinal and transverse velocity, 
longitudinal and transverse turbulent stress, and Reynolds shear stress at three 
locations in the wake. 

It is important to mention that, in the previous experiment, no hot-wire data was 
obtained in the region of reverse flow between x/C=1.0 and x/C= 1.25. Hot-wire 
anemometry is virtually unusable in this region due to the fact that hot-wires rectify the 
signals, which include the instantaneous forward velocity as well as the reverse velocity. 
Consequently, velocities measured by hot-wires in reverse flow regions appear higher 
than they actually are. 

Except for minor grid alignment errors in the y direction, most of the curves are 
in very good agreement. As expected however, the longitudinal turbulence stress 
curves at x/C= 1.25 (Figure 4.19) and the transverse turbulence stress curves at 
xC=1.25 (Figure 4.20) show considerable disagreement. In the longitudinal 
turbulence case, the hot-wire data shows more turbulence in the middle of the wake 
while in the transverse turbulence case, the hot-wire data shows less turbulence in the 
middle of the wake. This is again due to the fact that hot-wire measurements give a 
misleading velocity profile close to the trailing edge because of the reverse flow in that 
region [Ref. l:p. 70]. 

The hot-wire anemometer is quite accurate in regions where there is no reverse 
flow. However, in order to obtain reliable data on the mean velocities and turbulent 
stresses in the reverse flow regions and the near wake of this airfoil-spoiler system. 
Laser Doppler Velocimetry yields better results. 
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Figure 4.17 Comparison of LDV and hot-wire data for 
longitudinal velocity. 
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Figure 4.18 Comparison of LDV and hot-wire data for 
transverse velocity. 
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Figure 4.19 Comparison of LDV and hot-wire data for 
longitudinal turbulence stress. 
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Figure 4.20 Comparison of LDV and hot-wire data for 
transverse turbulence stress. 
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Figure 4.21 Comparison of LDV and hot-wire data for 
Reynolds shear stress. 
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E. ERROR ANALYSIS 



In this experiment, 1000 velocity samples were taken by the LDV in each 
direction (U and V) at every measurement location. The purpose in taking so many 
samples was to provide enough values in the fluctuating regions to ensure that an 
accurate statistical mean could be computed. Consequently, in the data reduction 
process all 1000 values were averaged to obtain values for U ave and V ave at each 
location. The equations in Table II were then used to obtain Reynolds shear stress, 
longitudinal and transverse turbulent stress at these locations (still using all 1000 
samples). 

To determine its effect, a different data processing method was tested on the raw 
data from several runs. In this procedure, only the samples which fell within one 
standard deviation of the mean were considered (approximately 67% of the 1000 
samples) in computing L' ave , V , L' rms and V rm$ . Again the equations in Table II 
were employed, just as before. This new process had little effect on U aye but V , 
changed considerably (from -1.3 to -11.0 at one location tested). However, the V ave 
values did not always decrease (one location changed from 8.7 to 16.9). Therefore, the 
changes were not due to the removal of a bias error. The longitudinal and transverse 
turbulence stress values decreased in all cases. This was expected since the biggest 
contributors to these values are the velocity terms which differ greatly from the mean. 

This new data processing procedure did not improve the data quality by 
removing any biases (based on the fact that some values of L' ave and V ave increased 
while others decreased on using the new process). As a matter of fact, this process had 
a negative effect on the data quality in that it neglects velocity extremes which, in turn, 
lowers the apparent turbulence intensity of the flow field. For this reason, the original 
data processing procedure (consider all 1000 samples) was used in this experiment. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

A two-component LDV system was set up for the purpose of measuring mean 
and fluctuating quantities of an airfoil with a deflected spoiler. The results of this 
experiment have shown that this system provides an efficient and effective means of 
flow field measurement, including regions of reverse flow and large velocity 
fluctuations. Comparison with hot-wire data demonstrates the LDV's superiority over 
the hot-wire in measuring in the near wake region. 

The results also clearly show two recirculating regions in the wake of the airfoil. 
One vortex originates at the trailing edge of the airfoil itself while the other originates 
at the trailing edge of the spoiler. The turbulence region is concentrated in the near 
wake and dissipates rapidly in the downstream direction. 

B. RECOMMENDATIONS 

In order to establish a more complete data base, it is recommended that more 
LDV measurements be obtained at other than a = 0° and 5 = 30°. These results are 
important to the modeling of this flow field. Also, to remove the possibility of tunnel 
blockage errors, it is necessary to perform further studies in a large scale wind tunnel. 

When using an LDV to measure a turbulent flow field such as in this experiment, 
a large number of samples should be taken at each location. Furthermore, unless 
attempting to remove a velocity bias all of the samples taken should be considered in 
computing L’ ave , V ave , etc. 
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APPENDIX 

COMPUTER PROGRAMS 



PURPOSE 



(1) TO CONVERT A 24 HOUR CLOCK TO A 12 HOUR CLOCK. 

(2) TO ADD A.M. OR P.II. TO THE "ADJUSTED" TIME . 



10 

20 

30 



40 

50 



IF 

CO 

IF 

IF 



CO TO 



CO 

CO 



TO 

TO 



30 

10 

30 

30 



SUBROUTINE STDTM ( I DATE* I TIME, I VAL) 
DIMENSION I DATE ( 3 ) 

DIMENSION IT l TIE (3 ) 

CALL DATE ( I DATE , I ERROR) 

CALL TIME ( IT I ME, I ERROR) 

IF ( ITIME( 1 ) .GE. 13) GO TO 
( IT1ME( l ) .EQ. 12) 

TO 20 

( ITIME<2) .CT.OO) 

( ITIME( 3 ) .CT.OO) 

I VAL = "AM" 

CO TO 40 
I VAL = "PM" 

IF ( ITI ME( 1 ) .LT. 13) 

I T l ME ( 1 ) = ITI ME ( 1 ) 

CO TO 50 

IF ( ITI ME ( 1 ) . CE . 0 1 ) 

ITlME(l) = ITI ME ( 1 ) 

CONTINUE 
RETURN 
END 



GO 



GO 

+ 



TO 
1 2 

TO 

12 



40 



50 
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c ****************************** 
c ****************************** 



c ** 




** 


G ** 


SPOILER.FR 


** 


C ** 




** 


c ** 


MIKE FOREMAN 


** 


c ** 




** 



g ****************************** 
c * * * * ************ * * * * * * * * * * * ** * 
0 

C LASER DATA ACQUISITION PROGRAM 



C 

G RLDR SPOILER STDTM DI024 GEL IRQ 



C 

c 

c 



c 



c 



c 

c 



DIMENSION I F I LE( 4 ) , lDT(tO), VAR(IO) 

DIMENSION ITI ME(3 ) * 1 I>ATE<3) 

DIMENSION I RAY 1 ( 1 000 ) , UUY2(I000) 

DIMENSION HUNK 16) • IBIN2(I6), VELI ( 1 000 ) , VEL2 (1000) 
DIMENSION V ( 5 , 20 ) , X ( 30 ) ,Y<30) 



INITIALIZE 

DO 10 1=1,3 
1 DT ( I ) = 0 
VA1K I ) = 10 . 

10 CONTINUE 

ACCEPT "RUN NUMBER (ZERO TO QUIT): ", IRUN 
I F ( IRUN.LE.O) CO TO 909 

15 ACCEPT "NUMBER OF SAMPLES PER COMPONENT: * , NUM 
I F ( NUM. LE . 1 000 . AND . NUM.GE . 0 ) GO TO 20 
TYPE "NUMBER OF SAMPLES MUST BE BETWEEN I AND 1000" 
CO TO 15 

20 ACCEPT "MIXING FREQUENCY *U " , FREQ 1 
ACCEPT "MIXING FREQUENCY *2: " , FIIEQ 2 



CONI =16.65 
C0N2= 16.65 
YAH ( l )=C0N1 
VAR(2)=C0N2 

SUBDEVICE ADDRESSES 



I SUB 1 =46K 
I SUB2=47K 
C 

C TIME PROGRAM 
C 

CALL STDTM ( 1 DATE , ITI ME , IVAL) 

C 

DO 30 1=1, 3 
I DT ( I ) = 1 DATE( I ) 

1 DT ( I +3 ) = ITI HE< I ) 

30 CONTINUE 

I DT ( 7 ) = I VAL 
I DIM 3) = I RUN 
YAU( 9 ) = FBEQ1 
VAR( 10)= FBEQ2 
C 

35 TYPE "WILL THIS BE A : (I) WAKE PROFILE" 

ACCEPT " (2) BOUNDARY LAYER PROFILE " 

C 

I F < IPROF.F.Q. I .OR. I PROF. EG. 2) CO TO 40 
TYPE "MUST ENTER I OR 2" 



, I PROF 
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GO TO 35 

40 I F ( I PROF. EQ. 2) GO TO 41 

ACCEPT M X-LOCAT 1 ON FOR THIS RUN (XT/C): - t VAR(3) 

GO TO 42 

41 ACCEPT "X-LOCATION FORTIUS RIJN (X/C): ",VAR(3) 

42 ACCEPT " NUMBER OF Y-LOCATIOHS: " , NflAX 
IF( lPROr.EQ.2) GO TO 45 

ACCEPT “STANDARD WAKE PROFILE Y-LOCATIONS? (1=YES,0=N0) " , LL 
1 F ( LL . EQ . 1 ) GO TO 44 

TYPE" INPUT Y-LOCATIONS IN INCHES (LOWEST POINT FIRST): - 
DO 43 1=1, NNAX 
WRITE (10,49) I 
ACCEPT*’ u , Y( 1 ) 

43 CONTINUE 
GO TO 50 

44 Y ( 1 )=-I .0 
Y(2)=-0 .8 
Y ( 3 ) =-0 . 6 
Y( 4 ) =-0 . 46 
Y ( 5 ) =-0 . 34 
Y ( 6 ) = -0 . 2 
Y ( 7 ) =-0 . 1 
Y( 8 ) =0 . 0 

Y ( 9 ) =0 . 1 
Y( 10) = .2 
Y( 1 1 )= .3 
Y ( 12) = .4 
Y( 13) = .5 
Y( 14)= .6 
Y( 15)= .7 
Y( 16)= .0 
Y( 17) = .9 
Y ( 1 0 ) = 1 .0 
Y ( l 9 ) = 1 .2 
Y ( 20 ) = 1 . 4 
GO TO 50 

45 TYPE ‘'STANDARD BOUNDARY LAYER PROFILE Y-LOCATIONS? ** 

ACCEPT " ( 1 =YES , 0=N0): u , LK 

1F(LK.EQ. 1 ) GO TO 47 

TYPE" INPUT Y-LOCATIONS IN INCHES ( LOWEST POINT FIRST):" 

DO 46 1=1 , NMAX 
WRITE ( 10,49) 1 
ACCEPT " " , Y ( I ) 

46 CONTINUE 
CO TO 50 

47 Y ( 1 ) = . 1 
Y(2)= .2 
Y ( 3 ) = .3 
Y(4)= .4 
Y ( 5 ) = .5 
Y(6)= .6 
Y(7)= .7 
Y(8) = .8 
Y(9)= .9 
Y( 1 0 ) = 1 .0 
Y ( 1 1 ) = 1 . 1 
Y ( 1 2 ) = 1 .2 
Y ( 1 3 ) = 1 .3 
Y ( 1 4 ) = 1 .4 
Y ( 1 5 ) = 1 .5 

49 F ORM \T ( 6X t "Y(“,I2,"):") 

50 KHITEdO, 125) 

125 FORMAT (" " , "BEGIN (SCANI .FR) >>",Z) 

CALL FSWAP( "SCANI .SV" ) 
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ACCEPT M D0 YOU WANT RAW DATA STORED? ( 1 =YES , 0 = NO) - f I RAW 
I DT ( 9 ) = I RAW 
I DT ( 10) = I PROF 
C 

C BEGIN TAKING DATA 
C 

N = 0 

57 CALL DIOI I ( ISUB1 , IRAY1 , IRAY2,N1JM, IER) 

C 

C DECODE 
C 

AVE1 =0 . 

AVE2=0 . 

60 DO 90 1=1 , NUM 

DO 65 J=1 t 16 

IBIN1 ( J ) = 1 AND ( 1 RAY 1 (I ) ,2**(J-1 ) ) 

IBIN2(J) = IANI)( I IIAY2 ( I ) ,2**!J-1 )) 

IF ( 1BIN1 ( J) .NE.O) IB1N1 ( J )=1 
I F ( I B I N2 ( J ) . NE . 0 ) I B1 N2 ( J ) = 1 
65 CONTINUE 
I BN 1 = 0 
I BN2=0 

DO 70 J = 0 , 9 

I BN 1 = I BN 1 +1 BI N 1 ( J+ 1 )*2**J 
I BN2= I BN2+ I B I N2 ( J+ 1 ) *2**J 
70 CONTINUE 

I EXP 1 =0 
I EXP2=0 

DO 75 J= 1 0 ♦ 13 

1 EXP 1 = 1 EXP 1 + I B1 N 1 (J + l ) *2** ( J- 1 0 ) 

1 EXP2= I EXP2+ I B I N2 ! J + 1 )*2**( J-10) 

75 CONTINUE 

C 
C 

VEL1 ( I ) = CON 1 * ( 2 . ** ( 5-1 EXP 1 )* 1 000 . /FLOAT! IBN1 )-40.+FREQl ) 

VEL2 ( I )=CON2*(2.**(5-IEXP2)*1000./FLOAT( I BN2 ) -40 . +FREQ2 ) 

AVE1 =AVE1 +VEL 1 ( I ) 

A V E2 = A V E2 + V EL2 ( I ) 

C 

90 CONTINUE 
C 

AVE1=AVE1/NUM 
AVE2= AVE2/NUM 
C 

C CHECK VELOCITIES ARE IN LIMITS 
C 

NN=N+1 

WRITFMO 92 1NN AVF2 AVF 1 

92 FORMAT! IX, "LOCATION: " ,12, 4X , M U ( AVE ) = " f F9 . 5 ♦ 4X , " V ( AVE ) = “,F9.5) 
1 F ! AVE1 . GT . 35 . ) GO TO 95 
IF!AVE2.GT. 140. ) GO TO 95 
GO TO 90 

95 WRITE FREE! 10) M BAD DATA: LOCATION " , NN 

ACCEPT " RETRY OR SKIP TO NEXT LOCATION? ! 1 =RETRY, 0=SKIP): " , ITRY 
I F ! I TRY . EQ . 0 ) GO TO 90 
CALL DIORM 
GO TO 57 
C 

90 IF! IRAW.NE. 1 ) GO TO 100 
C 

C STORE RAW DATA 
C 

CALL FOPEN ! 5 , "COUNTER" ) 

READ BINARY! 5) NJ 
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C 



C 



G 

C 



CALL FCLOSE( 5 ) 

NJ = NJ+ 1 

CALL FOPEN ( 5 , "COUNTER" ) 

WHITE BINARY ( 5 ) NJ 
CALL FCLOSE( 5 ) 

NF I LE = NJ 

IF(N.EQ.O) JFI LE=NF I LE 
I R=MOU ( NF I LE ,10) 

1 S=NF I LE/1 00 

IL= (NFI LE-IS*1GG)/I0 

I F I LE ( I ) = * RA * 

I F I LE( 2 ) = ’ VD’ 

I F I LE ( 3 ) = ( I S+60K) *400K+ ( IL+60K) 
I FI LE(4 ) = ( IR+60K)*400K 



CALL FOPEN (5 , I FILE) 
WRITE BINARY(5) (IDT(I) 
WRITE BI NARY ( 5 ) 

WRITE BINARY ( 5 ) 

WRITE BINARY ( 5 ) 

CALL FCLOSE ( 5 ) 



1 = 1 , 10 ) 

( VAR ( I ) , I = I , 10) ,NUM,Y(NN) 
( I RAY 1 ( I ) , I = 1 ,NUM) 

( 1 11AY2 ( 1 ) , 1 = 1 , NUM ) 



100 WRITE ( 10,99) 

99 FORMAT <" ","*** BEGIN RMS CALCULATIONS >> 
URMS=0 . 

VIUIS = G. 

UVRMS=G . 

DO 110 1=1 , NUM 

URTIS = URMS+( ( VEL2 ( I ) )**2- ( AVE2 ) **2 ) 

vims=viuis+( (veli < i )**2 -<avei )**2) ) 

UVRIiS = UVRTIS + ( VEL 1 ( I )^VEL2( I ) ) - ( AVE1 *AVE2 ) 
110 CONTINUE 

URIiS = (UimS/FLOAT(NUN) ) 

VRMS= ( VRTIS/ FLOAT ( NUN ) ) 

UVRI1;5= ( UVRTIS/FLOAT ( NUT1 ) ) 

WRITE (10,111) 

I 1 1 F0RI1AT ( " + " f u > > END RTIS CALCULAT I ONS " ) 

N = n+ 1 

V( 1 , N ) = AVE2 

V ( 2 , N ) = A V E 1 
V ( 3 , N ) =URHS 

V ( 4 , N ) = V Rf IS 

V ( 5 , N ) =UVIU1S 

I F ( N . CE . NI'LiX) GO TO 120 



MOVE TRAVERSE 



Z) 



CALL FSWAP( "TRELAY.SV" ) 

ACCEPT "TRAVE1RSE MOVED... ANY NUMBER TO CONTINUE. " 



C 

CALL DIORII 
CO TO 57 
C 

C OBTAIN TEMP, STAG PRES., MACH *, ETC 



120 I F ( IRUN.EQ. 1C00) GO TO 127 
CALL FSWAP( "TRELAY.SV" ) 

WRITE ( 10,126) 

126 FORMAT ("TRAVERSE RESET TO INITIAL POSITION") 

127 CALL FOPEN (5, "SCAN I " ) 

READ Bl NARY ( 5 ) VAR ( 4 ) , VAR ( 5 ) , VAR ( 6 ) , VAR ( 7 ) , VAR( 8 ) 
CALL FCLOSEC 5 ) 



C 

C WRITE DATA INTO PERMANENT DATA FILE 



, IAA 
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CALL FOPEN ( 5 , "PERMCOUNT" ) 

READ Li I NARY ( 5 ) NF 
CALL FCLOSE ( 5 ) 

WRITE FREE(IO) " LAST DATA FILE USED: " ,MF 
WHITE FREE(IO) "RUN NUMBER: " , 1RUN 

ACCEPT "’DATA* FILE NUMBER (NEC * TO SKIP): " , MF I LE 
IF(MFILE.LT.O) CO TO 599 
CALL FOPEN ( 5 , " PLIUICOIJNT" ) 

WRITE BINARY ( 5 ) NFILE 
CALL FCLOSE( 5 ) 

I R = MOD( FIF I LE ,10) 

IL=IIFI LE/10 
I F I LE ( 1 ) = ’ DA * 

I F I LE( 2 ) = ’ TA * 

IFILE(3) = ( IL+60IC)*400IC+( IR+60K) 

I F I LE ( 4 ) = 0 

CALL FOPEN (5, I FILE) 

WRITE B I NARY ( 5 ) ( I DT (I ), 1 = 1 ,10) 

WRITE BINARY (5) (VAIK I ) , I = 1 , 1 0 ) , NUM , NMAX , JF I LE , NF I LE 
WRITE BI NARY ( 5 ) ( (VC I ,J) ,J=1 , NMAX) ,1 = 1,5) 

WRITE BI NARY ( 5 ) ( Y( 1 ) , I = I , NMAX) 

CALL INCLOSE ( 5 ) 

599 WRITE FREE ( i 0 ) "AMBIENT TEMP: 

WRITE FREE(IO) "ATMOSPHERIC PRESSURE: 

WRITE FREE (10) "STAGNATION PRESSURE: 

WRITE FREER 10) "FREE STREAM MACII NUMBER: 

WRITE FREE (10) "FREE STREAM VELOCITY: 

DO 650 1=1, NMAX 
Y ( I ) =Y( 1 )/6. 

650 CONTINUE 

ACCEPT" DO YOU DESIRE A CURVE? ( 1 =YES , 0 = N0) : " , I CURVE 
I F ( I CURVE . EQ . 0 ) CO TO 973 
WRITE FREE(IO) 

600 TYPE "WHAT PROFILE DO YOU DESIRE?" 

TYPE" (1) U(AVE)/U( INFINITY)" 

TYPE" (2) V(AVE)/U( INFINITY) " 

TYPE" (3) U ( RMS )/ ( U ( I NF I N I TY) **2 ) " 

TYPE" (4) V(RMS)/(U( I NFI N 1TY>**2 ) " 

ACCEPT" (3) UV(1UIS)/(U( INFINITY)**2) ",ISEL 

B I GX=0 . 

DO 700 1 = 1, NTIAX 

I F ( ISEL.LE.2) X( I )=V( I SEE, I )/VAR(0) 

I F ( ISEL.CT.2) X( I ) = V ( ISEL, I )/ ( VAR( 3)**2 ) 

IF(X( I ) .CT.BICX) BI CX=X( I ) 

700 CONTINUE 

CALL FOPEN (5, "PROFILE" ) 

WRITE BI NARY ( 5 ) NMAX, I PROF , I SEL , BI CX 
WRITE BI NARY ( 5 ) (X( I ) , I = 1 , NMAX) 

WRITE B I NARY ( 5 ) ( Y ( I ) , I = I , NMAX) 

CALL FCL0SE(5) 

CALL r SWAP ( " PLOTP . SV " ) 

ACCEPT" ANY * TO CONTINUE " ,ANY 
TYPE "DO YOU DESIRE ANOTHER CRAPII? " 

ACCEPT" ( 1 = YES , 0=N0 ) : " , I GR 

IF( ICR.EQ. 1 )C0 TO 600 
975 WRITE FREE (10)" ENSURE PRINTER IS ON" 

WRITE ( 12,980) I DATE , IT I ME, I VAL 

900 FORMAT ( 1 OX, "DATE: " ,2(12,"/" ) , 14, 5X, "TIME: " ,2(12," : " ) , 12, IX, A2) 
WRITE FREE (12) 

WRITE FREE( 1 2 ) "RUN NUMBER: ",IRUN 
WRITE FREE( 1 2 ) 

1 F ( I PROF .EQ. 1 ) WRITE FREE( 1 2 ) " WAKE PROFILE" 



" , VAR( 4 ) , " DECF" 
" , VAR ( 5 ) , " PS I" 

" , VAR ( 6 ) , " PSI " 

" , VAR ( 7 ) 

" , VAR( 8) , " M/S" 
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I F ( irnOF.EQ.2) WRITE FREE (12)" BOUNDARY LAYER PROFILE" 

WRITE FREE (12)" X-LOCATION: ",VA1U3) 

WRITE FREEU2) 

WRITE FIIEE( 12) "NUMBER OF Y-LOCATIONS TAKEN: " , NMAX 
WRITE FREE( 1 2 ) 

WRITE( 12,981 ) 

981 F0RMAT(2X, "LOCATION" , 10X, "U(AVE) " , IOX, "V(AVE) " , 10X, "U(RMS) " , 

81 IX, "V(RTIS) " , 1 IX, "UV(RNS) " ) 

DO 935 1=1, NMAX 

WR I TE (12, 984 ) Y ( I ) , V( 1 , I ) , V(2 , I ) ,V<3 , 1 ) , V( 4 , I ) ,V( 5 , I ) 

984 FORMAT ( 2X , F8 . 6 , 9X, F9 . 4 , 8X , F9 . 4 , 7X, FI 0 . 4 , 7X , F 1 0 . 4 , 7X , F 1 O . 4 ) 

983 CONTINUE 

WRITE FREE( 12) 

WRITE FREE( 12) "NUMBER OF SAMPLES PER Y-LOCATION: ",NUM 
WRITE FREE( 12) "MIX FREQ 1 = " , FREQ I , " MIX FREQ **2= " , FREQ2 
WRITE FREE( 1 2 ) " VEL CONST *1= ",CONI," VEL CONST *2 = ",C0N2 
WRITE FREE( 1 2 ) 

I F ( I RAW . EQ . 0 ) WRITE FREE(12)"RAW DATA HAS NOT BEEN STORED" 

IF( IRAW.EQ. 1 ) WRITE FREE( 1 2 ) “RAW DATA IS STORED IN ’ RAWD* FILES ", 
8JFILE," THRU " , NF1 LE 

WRITE FREE( 12) "PERMANENT DATA FILE: DATA " * MF I LE 
WRITE FREE( 1 2 ) 

WRITE FREE( 12) "AMBIENT TEMP= ",VAR(4), n DECF " 

WRITE FREE( 12) "ATMOSPHERIC PRESSURE= ",VAR(5),“ PS I" 

WRITE FREE( 12) "STAGNATION PRESSURE^ M ,VAR(6)," PSI " 

WRITE FREE( 12) "FREE STREAM MAC11 * = " ,VAR(7) 

WHITE FREE( 12) "FREE STREAM VEL0C1TY= M ,VAR(0)," M/S" 

999 STOP 
END 
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xxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXXXXXXXX 



C xx 




XX 


C xx 


" SCAN I .FR" 


XX 


C xx 




XX 


C xx 


MIKE FOREMAN 


XX 


C xx 




XX 



xxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXXXXXXXX 

"SCAN I . Fll" TAKES T1IE DATA FROM THE 6 

SCAN 1 VALVES AND WRITES THIS DATA (IN COUNTS) 

ON THE SCREEN AND INTO A DATA FILE. 

DIMENSION I RAY ( 1 200 ) , Y<6,0:23) 

THE SCAN l VALVES ARE CONNECTED TO CHANNELS 1 THROUGH 6: 

N START = 1 
LIMIT = 6 
NUM = L 1 M I T-NSTART + 1 

ZEROIZE THE MATRIX 

DO 10 I = 1*6 
DO 10 J = 0,23 
Y( I , J ) = 0 
10 CONTINUE 

SCAN I VALVE RELAY INFORMATION 
I DEV ICE = DEVICE NUNDEU (25IC FOR SCAN I VALVE) 

RELAY NUMBERS: HOME = RELAY 

STEP = RELAY *2 

TO CONVERT REIAY NUMBER TO BINARY WORD USE: 

N=2**( 16-N) 



I DEV ICE = 25K 

N 1 =2 ; STEP = RE LAY-* 2 

ISTEP = 2** ( 1 6-N 1 ) 

N2 = 3 ; 1I0I1E = RELAY* 0 

I HOME = 2**(16-N2) 

RESET ALL BITS TO ZERO 



CALL DOUTW ( I DEV I CE , 0 ) 

SETUP THE A/D DEVICE WITH THE STARTING AND ENDING CHANNELS 



CALL FADST( NSTART , LIMIT , IER) 

RESET THE SCAN I VALVE TO THE HOME POSITION 
(ENERGIZE THE RESPECTIVE RELAY, THEN DEENERGIZE IT) 



C 

C 

c 



CALL DOUTW ( I DEV ICE, inOME) 

CALL FDELY ( 5 ) 

CALL DOUTW ( I DEV I CE , 0 ) 

TAKE 200 SAMPLES AT EACH PORT 

DO 30 L=0 , 23 
DO 20 N = 1 ,200 

CALL FADC'K I RAY ( ( N- 1 )*NUM+1 ), NUM, IER) 
20 CONTINUE 
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STEP THE SCAN I VALVES TO THE NEXT PORT 

CALL DOUTV( I DEV ICE, I STEP) 

CALL FDELY ( 5 ) 

CALL DOUTW ( I DEV I CE , 0 ) 

COMTUTE THE TOTAL COUNTS FOR EACH PORT (OF EACH SCAN I VALVE) 

AND SIORE EACH TOTAL IN THE ARRAY , Y. 

DO 30 M=NSTART , LIMIT 
Ml =M-NSTART+ 1 
DO 30 N= 1,206 

Y ( Ml , L)=Y(Ml ,L)+IRAY( (N-l ) *NUM+ ( M-NSTART+ 1 ) ) 

30 CONTINUE 

RETURN Til E SCAN I VALVES TO THE HOME POSITION 

CALL D0UTW( I DEV ICE, I HOME) 

CALL FDELY(5) 

CALL DOUTW ( I DEVICE, 0) 

D I SPLAY THE AVERAGE SCAN I VALVE VALUES (IN COUNTS) ON THE SCREEN 
NOTE: THE AVERAGE IS COMPUTED INSIDE THE WRITE STATEMENT ITSELF 



WRITE FREE ( 1 0 ) M BEGIN TEMP . SV . . . " 

CALL FSVAP ( “TEIIP . SV K ) 

WHITE FREE( I 0 ) “ END TEIIP.SV >>> BEGIN CALPRES . SV . . . " 

CALL ESWAP( “CALPRES. SV“ ) 

WHITE FREE ( 1 0 ) M END CALPRES. SV >>> BEGIN NPARO.SV... 

CALL FSWAP( ■NPARO.SV" ) 

WHITE FREE ( I 0 ) “ END NPARO.SV >>> CONTINUE SCANI.SV.. 

CALL T OPEN ( 5 , "TEMP “ ) 

READ BI NARY ( 5 ) TEMP 
CALL FCLOSE( 5 ) 

CALL FOPEN (5, -CALPRES" ) 

READ B I NARY ( 5 ) CAL 
CALL TCLOSE(5) 

CALL I OPEN ( 3 , "ATMPRES n ) 

READ BINARY (3) ATM 
CALL ECLOSE(3) 

ATfI=ATf 1^3377/639 3 ; CONVERT TO PSI 

WRITE EREE( 10) "CAL" , CAL 

TEfiPR = TEnP-r4 59.6 ; CONVERT TO RANK I NE 

TuTPHES= (Y( 1,2) -Y (1,0) )*CAL/(Y( 1 , 1 )-Y( 1 ,0) ) +ATM 
STATPRES = ( Y ( 4 , 2 ) - Y ( 4 , 0 ) K CAL/ ( Y ( 4 , I ) - Y ( 4 , 0 ) ) +ATM 
F:iAClI = S0RT(ABS(3.*( ( ( STATPRES/TOTPRES ) ** ( -2 . /7 . ) )-l . ) ) ) 

TRY1 = 5 . *( ( ( STATPRES/TOTPRES ) ** ( -2 . /7 . ) )-I . ) 

TRY2=ABS(TRYI ) 

IF(TRY1 .EQ.TRY2) GO TO 40 

WRITE FREE( 10) “STATPRES IS GREATER THAN TOTPRES M" 

40 T l =TEMPR/ ( 1 + ( FMACH**2 . ) /o . ) 

A = SORT ( I .4*171 6. *Tl )*.3040 

FVEL=A*FMACH 

CALL FOPEN ( 5 , "SCANI " ) 

WRITE BI NARY ( 5 ) TEMP, ATM, TOTPRES, FMACH , FYEL 
CALL FCLOSE ( 5 ) 

C REMOVE THE DEVICES 
C 

CALL DIORM 
CALL FADRM 
C 

SIOP 

END 
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c *************************** 
c *************************** 



c ** ** 

c ** "TEMP. Fir ** 

c ** ** 

c ** MiKE FOREMAN ** 

G ** ** 



G *************************** 
C *************************** 



G 

C "TEMP . FR" COMPUTES TIIE CURRENT ATMOSPHERIC 
C TEMPERATURE AT TIIE INLET TO TIIE INDRAFT 
G TUNNEL . THIS VALUE IS OUTPUT TO THE DISPLAY 
C SCREEN AND ALSO WRITTEN INTO A FILE CALLED 
C "TEMP". 

C 

DIMENSION I RAY ( 200 ) 

C 

C TIIE TEMPERATURE PROBE INTERFACES WITH CHANNEL 8 
G 

N START = 8 
LIMIT = 8 

NUM = L I M I T-NSTART+ I 
C 

C SET UP THE DEVICE 
C 

CALL FADST(NSTART, LIMIT, IER) 

C 

C TAKE 200 SAMPLES 
C 

DO 10 1=1 , 200 

CALL FADCT ( I RAY ( I ) , NUM , I ER) 

10 CONTINUE 
C 

C COMPUTE TIIE TOTAL NUMBER OF COUNTS FOR THE 200 SAMPLES 
G 

TEHP=0 . 

DO 20 1=1 ,200 
TEMP = TEMP + 1 RAY ( I ) 

20 CONTINUE 
C 

G CALCULATE TIIE AVERAGE OF TIIE 200 SAMPLES 
C 

TEIIP=TEMP/200. 

G 

C THE TEMPERATURE SIGNAL ENTERS THE AMPLIFIER IN VOLTS. 

G TIIE AMPLIFIER GIVES TIIE SIGNAL A 60DB GAIN (FACTOR 

C OF 1000). THEREFORE , TIIE SIGNAL GOES TO TIIE A/D CONVERTER WITH DIMENSIONS 
G OF MILLIVOLTS. AFTER CONVERSION, TIIE FACTOR 1638.3 COUNTS/MILLIVOLT 
C IS NECESSARY TO CONVERT THE COMPUTER UNITS (COUNTS) BACK TO MILLIVOLTS: 

C 



TEIIP =TEIIP/ 1 638 . 3 ; CONVERT COUNTS BACK TO MILLIVOLTS 

C 

C TIIE "TEMP" HAS NOW BEEN CONVERTED TO MILLIVOLTS 
C HOWEVER, TIIE TEMPERATURE EQUATION REQUIRES THAT 
C THE "TEMP" BE IN MICROVOLTS: 

C 

TEMP=TEMP* 1 000 . 

C 

C 

C BECKMAN COEFFICIENTS 
G 

CO = 3 . 19227G7E+01 
Cl =3 . 5497620E-02 
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nno onnn ono 



C2= -2 . 4786653E-07 



CALCULATE TEMPERATURE 

TEMP=C0+C1*TEMP+C2*TEMP**2 
WRITE TEMP TO SCREEN AND FILE "TEMP M 
WRITE FREE (10) 

WRITE FREE (10) "TEMPERATURE ( DEGF ) = M ,TEMP 
CALL FOPEN (5, "TEMP" ) 

WRITE BINARY (5) TEMP 
CALL FCLOSE (5) 

DISCONNECT THE DEVICE 

CALL FADRM 

STOP 

END 
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c *************************** 
c *************************** 



G ** 




** 


G ** 


" CALPRES. FR" 


** 


c ** 




** 


C ** 


MIKE FOREMAN 


** 


G ** 




** 



c *************************** 



o *************************** 

c 

C "CALPRES . FR" THIS PROGRAM COMPUTES THE 
C CALI BRATEI) PRESSURE TO TIIE SCAN I VALVES. 

C 

DIMENSION I RAY ( 200 ) 

C 

C THE CALIBRATION PRESSURE IS CONNECTED TO CHANNEL 7. 
G 

N START = 7 
LIMIT = 7 

Nim = L I N I T-NSTART+ 1 
C 

C SET UP TIIE A/D DEVICE: 

0 

CALL FADST( NSTART , LIMIT , IER) 

C 

C TAKE 200 SAMPLES AT CHANNEL 7 



DO 10 1=1 ,200 

CALL FADCT( IRAY( I ) , NUM, IER) 

10 CONTINUE 
G 

G COMPUTE THE AVERAGE OF TIIE 200 SAMPLES 
C 

Y = 0 . 

DO 20 1=1 ,200 
Y=Y+FLOAT ( I RAY ( I ) ) 

20 CONTINUE 
Y= Y/2C0 . 



C 

G CALIBRATION PRESSURE IN INCHES OF HG = VOLTS* 10 



CAL= < Y/ 1 638 . 3 ) * 1 0 . 

G 

G CALI BILAT I ON PRESSURE IN PS I 



INCHES OF HC * (3377/6895) 



C 

CAL=CAL* ( 3377 . /6895 . ) 

G WRITE FREE(IO) CAL 

C 

G WRITE THE CALIBRATION PRESSURE TO THE SCREEN AND TO 
C A TEMPORARY FILE CALLED "CALPRES", 

C 

G WRITE FREE(IO) 

G WRITE (10,25) CAL 

25 F ORMAT ( " CAL I BRAT I ON PRESSURE: ",F7.4, M PSI") 

G 

CALL FOPEN ( 5 , "CALPRES " ) 

WRITE BINARY (5) CAL 
CALL FCLOSE (5) 

G 

G REMOVE THE DEVICE 
C 

CALL FADRM 

STOP 

END 
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C “ PLOTP . FR" IS CALLED BY ”SPO I LER. FR n TO PLOT SELECTED VELOCITY 
C PROFILES. 

C 

DIMENSION X( 25 ) , Y(25) 

DIMENSION IST1 (3) , IST2(8) , IST3(8) , IST4(8) 

D I MENS I ON I ST 5 ( 8 ) , I ST 6 ( 9 ) , I ST7 ( 2 1 ) , I ST8 (31) 

C 

CALL CURS I Z( 4) 

DATA I ST 1/09, 47 ,67/ 

DATA I ST2/05 , 40 ,65,06,69,41 ,47,85/ 

DATA I ST3/86 ,40,65,86,69,41,47, 05/ 

DATA I ST4/C5 ,40,82,77,83,41,47, 85/ 

DATA I ST5/06 ,40,02,77,03,41,47, 05/ 

DATA I ST6/85 , 86 , 40 , 82 , 77 , 83 , 41 ,47,05/ 

DATA I ST7/C7 ,65,75,69,32,06,69,76,79,67,73,84,89,32, 
*80,02,79,70,73,76,69/ 

DATA I ST8/66 ,79,05,78,60,65,02,09,32,76,65,89,69,82,32, 
*06,69,76,79,67,73,04,09,32,80,82,79,70,73,76,69/ 

C 

CALL F0PEN(5, “PROFILE" ) 

READ BINARYC5) NMAX, I PROF , ISEL, BIGX 
READ B I NARY ( 5 ) ( X( 1 ) , I = 1 , NILAX) 

HEAD BI NAIIY ( 5 ) (Y(I), 1 = 1, NIlAX) 

CALL FCLOSE ( 5 ) 

I F ( B l GX . LT . .5) BI GX= . 5 
MUTE FREE ( I 0 ) Y( 1 ) ,Y( NIlAX) ,BIGX 
ACCEPT "ANY NUIIBLR TO CONTINUE: " , IRAND 
C 

CALL 1 H I TT ( 960 ) 

CALL B I N I TT 
CALL CililS I Z( 3 ) 

CALL TERM ( 3 , 1023) 

CALL SLI MX (ICO, 900 ) 

CALL SLHIY( 150,650) 

CALL XFiU I ( 5 ) 

CALL YFIUK5) 

CALL XTICS( 10) 

CALL YTICS( 10) 

CALL DLI MX( 0 . 0 , BIGX) 

CALL DLI MY ( Y( 1 ) , Y ( NMAX) ) 

CALL NP I S (NILAX) 

CALL C11ECE(X,Y) 

CALL DSPLAY (X , Y) 

CALL C11HS I Z ( 3 ) 

CALL IiOVALS( 10,525) 

CALL VLABEL (3, IST1 ) 

CALL MOVABS (300,725) 

IF( IPROF.EQ. 1 ) CALL HLABEL(2I , IST7 ) 

I F ( I PROF . EQ . 2 ) CALL HLABEL(3I , IST8) 

CALL MOVABS (420, 75) 

IF( ISEL.EQ. 1 ) CALL HLABEL(8 , IST2) 

I F ( ISEL.EQ.2) CALL HLABEL( 8,1 ST3) 

I F ( ISEL.EQ. 3) CALL HLABEL (8,1 ST4 ) 

I F( ISEL.EQ. 4) CALL IILABEL ( 0 , I ST5 ) 

I F ( ISEL.EQ. 5) CALL IILABEL ( 9 , I ST6 ) 

CALL F I N I TT (0,7 50 ) 

CALL C11RS 1 Z ( 4 ) 

S'FOP 

END 
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******************************* 

******************************* 



C ** 




XX 


c ** 


II I STOGRAM . FR 


** 


c ** 




XX 


c ** 


MIKE FOREMAN 


XX 


c ** 




XX 



****************** ************* 
******************************* 



iiistograii.fr records data from the ldv 

COUNTER AND PLOTS A HISTOGRAM OF THIS DATA. 

REQUIRED SUBROUTINES: DI022 uFLIBG 

DIMENSION I RAY ( 1000) , IDIN( 16) , FREQ( 1 000 ) 

DIMENSION X(101),Y(101) ,VKL< 1000) 

TYPE “INPUT SUBDEVICE CODE TWO DIGITS ONLY : " 

TYPE " CHANNEL 1 = 46“ 

TYPE " CHANNEL 2 = 47" 

READ( 11,5) ISUB 
5 FORMAT (02) 

ACCEPT “VELOCITY CONSTANT? " ,CON 
ACCEPT "MIXING FREQUENCY? " , FREQM 
10 ACCEPT " NUMBER OF SAMPLES (0 TO STOP) ?",NUM 
IF (MJM.GT. lOOO.OR.NUfl.LT.O) GOTO 10 
IF (NUri.EQ.O) CO TO 990 
CALL 1)1 01 I (ISUB, l RAY, NUM, I EH) 

DO 100 J=1 , NUM 
DO 50 1=1,16 

I BI N ( I ) = I AND( I RAY ( J ) ,2**( 1-1 ) ) 

50 I F ( I B I H ( I ) . N E . 0 ) I B I H ( I ) = 1 

I BN = 0 

DO 60 1=0,9 
I BN= I BN+ 1 B I N ( 1+1 )*2**I 
60 CONTINUE 

I EXP = 0 

DO 70 1=10,13 

1 CXP= I EXP + 1 B I N ( I + i )*2**( I -10) 

70 CONTINUE 

FREQ ( J )= 32000 . /FLOAT ( I BN*2** I EX? ) 

VEL( J ) =CON* ( FREQ ( J ) -40 . + FREQM ) 

100 CONTINUE 

AVER=0. 

A V E = 0 . 

DO I 1 0 1 = 1, NUM 
AVEll=AVER+FREQ( I ) 

AVE=AVE+VEL( I ) 

I 10 CONTINUE 

AVER= AVER/NUM 
AVE= AVE/NUM 

WRITE FREE(IO) "AVERAGE DOPPLER FREQUENCY= ", AVER, “MHZ" 
WRITE FIIEE(IO) "AVERAGE VELOCITY* ",AVE," M/S" 

ACCEPT “DO YOU WISH A H I STOGILAM? (YES=1, N0 = 0 ) " , AA 
IF(AA.EQ.O) GO TO 999 
NM= I F I X ( AVE ) 

J J = MM-5 

lF(JJ.LT.O) JJ=0 
KK = J.I + 1 0 
LL= UJ*10)-1 
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IV=0 

z=o. 

DO 200 1=1,101 
K= I +LL 

X( I ) = ( FLOAT (K) )/I0. 

y ( i ) =o . 

DO 150 J=1 , NUM 
I 1 = 1 F I X ( 10.*VEL( J) ) 

I F ( 1 I .EQ.K) V( I ) =Y( I ) + I 
150 CONTINUE 

WRITE FREE (10) X( I ) ,Y( I ) 

I W= 1 V+ 1 F 1 X ( Y ( I ) ) 

1 F ( Y( I ) . GT . Z ) Z= Y ( 1 ) 

200 CONTINUE 

Z= FLOAT ( ( 1F1XCZ/100. )*10O)+1O0) 

IFCZ.CT. 1000. ) Z=1000. 

AB= FLOAT ( JJ) 

AC = FLOAT (KK) 

WRITE FREE (10) AB,AC,Z 

CALL FOPEN ( 5 , **H 1 STO" ) 

WRITE BINARY (5) AB,AC,Z 
WRITE BINARY (5) (X(l), 1=1,101) 

WRITE BINARY (5) (Y(I), 1=1,101) 

CALL F CLOSE (5) 

CALL FSWAP ( "HISTOGRAPH.SV" ) 

WRITE FREE(IO) 

WRITE FREE( 10) 

WRITE FUEE(IO) 

WRITE FREE(IO) 

WRITE FREE (10) " TOTAL SAMPLES PLOTTED: IV 

CALL DIORM 
999 STOP 
END 
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C MflSTOGRAPH.FR" USED BV II I STOCRAM . FR TO PLOT 
C A 11 1 STOGllAM OF LASER DATA. 

C 

DIMENSION X( 101 ) ,Y( 101 ) 

DIMENSION 1 ST 1(7), 1ST2(21), IST3(15) 

C 

DATA 1 ST 1 /G3 ,65,77,00,76,69,03/ 

DATA 1 ST2/72 ,73,03,04,79,71,82,65,77, 

*32 ,79,70,32,76,60,06,32,60,65,84,65/ 

DATA 1 ST3/7G , 82 ,69,01 ,85,69,70,67,89,32,40,77,72,90,41/ 
C 

CALL FOPEN ( 5 , "HI STO '* ) 

READ BINARY (5) AB,AC,Z 
READ BINARY (5) (X<1),I=I ,101) 

READ BINARY (5) ( Y ( 1 ) , I = 1 , 1 0 1 ) 

CALL FCLOSE (5) 

C 

G 

CALL 1 N ITT ( 960 ) 

CALL BlNITf 
CALL TERM (3, 1023) 

CALL VRARST( 15.6,3) 

CALL S L 1 MX ( 100,900) 

CALL SL1I1Y ( 150,650) 

CALL XFRM ( 2 ) 

CALL YFIIH ( 5 ) 

CALL XTICS< 10) 

call yricsc to) 

CALL XfiTCS (10) 

CALL YMTCS(2) 

CALL DLI MX( AB , AC ) 

CALL DLI MY ( 0 .0 , Z) 

CALL NETS (101) 

CALL CHECK (X, Y) 

CALL DSPLAY ( X , Y ) 

CALL MOV A BS< 10,525) 

CALL VLABEL(7 , ISTl ) 

CALL MOVABS (350,725) 

CALL I1LABEL( 2 1 , I ST2 ) 

CALL MOVABS (420, 75) 

CALL 11 LA BEL ( 15, 1 ST3 ) 

CALL MOVABS (900, 150) 

CALL BIUVABS ( 900 , 650 ) 

CALL CURS 1 Z( 4 ) 

CALL F 1 N 1 TT( 0 , 750 ) 

STOP 

END 
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